Strong magnetic coupling of an inhomogeneous NV ensemble to a cavity 
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We study experimentally and theoretically a dense ensemble of negatively charged nitrogen- 
vacancy centers in diamond coupled to a high Q superconducting coplanar waveguide cavity mode 
at low temperature. The nitrogen-vacancy centers are modeled as effective spin one defects with 
inhomogeneous frequency distribution. For a large enough ensemble the effective magnetic coupling 
of the collective spin dominates the mode losses and inhomogeneous broadening of the ensemble 
and the system exhibits well resolved normal mode splitting in probe transmission spectra. We 
use several theoretical approaches to model the probe spectra and the number and frequency dis- 
tribution of the spins. This analysis reveals an only slowly temperature dependent q-Gaussian 
energy distribution of the defects with a yet unexplained decrease of effectively coupled spins at 
very low temperatures below 100 mK. Based on the system parameters we predict the possibility 
to implement an extremely stable maser by adding an external pump to the system. 



PACS numbers: 42.50.Pq,42.50.Ct,61.72.jn,03.67.-a 

I. INTRODUCTION 

Systems of spin ensembles coupled to a cavity mode 
are considered a promising physical realization for pro- 
cessing and storage of quantum information [I], ultra- 
sensitive high resolution magnetometers [2] or localized 
field probes. Collective magnetic coupling of a large 
ensemble to the field mode allows to reach the strong 
coupling regime, even if a single particle is hardly cou- 
pled. Implementations based on superconducting copla- 
nar waveguide (CPW) resonators provide an interface of 
the spin ensemble with processing units for supercon- 
ducting qubits. Here the spins can serve as a quantum 
memory or as a bridge to optical readout and communi- 
cation j^. 

Different types of ensembles were proposed for this 
setup, ranging from clouds of ultracold atoms [3] over 
polar molecules [1] to solid state systems like rare-earth 
spin ensembles [5] or color centers in diamond [6l [J • Here 
we focus on the negatively charged nitrogen-vacancy de- 
fects in diamond (NV). Those are naturally present in 
diamond but can also be readily engineered with very 
high densities, still maintaining long lifetimes and slow 
dephasing in particular at low temperatures of T < 1 K. 
In the optical domain, they are extremely stable and very 
well studied since many years [5]. 

The magnetic properties of the relevant defect states 
can be conveniently modeled by effective independent 
spin 1 particles, where the effective local interaction of 
the electrons within the defect shifts the (mg = ±1) 
states with respect to the ms — state On the 
one hand, this coupling provides the desired energy gap 
in the 3 GHz regime, but on the other hand, as a con- 
sequence of local variations of the crystal field, this shift 
exhibits a frequency distribution leading to an inhomoge- 
neous broadening of the ensemble. The inhomogeneities 
are thought to be predominantly caused by crystal strain 
and excess nitrogen, which is not paired with a neighbor- 



ing vacancy [TU] . 

While for a perfectly monochromatic ensemble of N 
particles in a cavity, strong coupling simply requires an 
effective coupling gsingieV^ larger than the cavity and 
spin decay rates, not only the width, but also the details 
of the inhomogeneous distribution are known to strongly 
influence the dynamic properties of the real world sys- 
tem [llH13j . In particular a Gaussian or a Lorentzian dis- 
tribution of equal half width, lead to different widths and 
magnitudes of the vacuum Rabi splitting. Only above a 
critical coupling strength the rephasing via common cou- 
pling to the cavity mode will prevent dephasing of the 
collective excitation and lead to a well resolved vacuum 
Rabi splitting. 

In our theoretical studies of this system we use differ- 
ent approximation levels to analyze the central physical 
effects present in an inhomogeneously broadened system, 
as they are observed in the measurements. While many 
qualitative features can be readily understood from a 
simple, coupled damped oscillator model with an effective 
linewidth, a detailed understanding of the observed fre- 
quency shifts and coupling strengths in the experiments 
requires more sophisticated modeling of the energy dis- 
tributions and dephasing mechanisms. In particular the 
observed temperature dependence relies on a finite tem- 
perature master equation treatment of a collective spin 
with proper dephasing terms. This is compared to the 
experimental values of up to {N w 10^^) particles as pre- 
sented already in [7j- 

The paper is organized as follows: The general prop- 
erties of the system are introduced in Sec. [TTj A first 
approximative treatment using coupled harmonic oscil- 
lators is shown in Sec. |III A} In Sec. |IIIB] we incorporate 
the inhomogeneous frequency distribution of the NVs via 
an extra decay of the polarization of the ensemble. In 
this context we also analyze the effects of thermal exci- 
tations in the spin ensemble and in the cavity. Finally, 
in Sec. |III C] we interpret our measurements using the re- 
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solvent formalism in order to extract the exact form of 
the inhomogeneity |13j . The prospect of implementing 
a narrow bandwidth transmission line micro-maser using 
an inhomogeneously broadened ensemble is discussed in 
SecHVl 



We thus approximate the composed system of cavity and 
ensemble with the Tavis-Cummings Hamiltonian 



II. GENERAL SYSTEM PROPERTIES AND 
EXPERIMENTAL IMPLEMENTATION 

The ground state of the NV center is a spin triplet 
(5 = 1), where the ms — ±1 states are split from 
the ms = state by about 2.88 GHz at zero mag- 
netic field [in]. In addition the degeneracy between the 
ms = ±1 states is lifted due to the broken C31, symmetry. 
Applying a homogeneous magnetic field enables Zeeman 
tuning of the TO5 — ±1 states, so that the ms = — > ±1 
transitions can be tuned selectively into resonance with 
the CPW cavity. 

Note that here the Zeeman tuning is also varying with 
the NV centers orientation relative to the applied field 
direction. In NV center ensembles all the four symmetry 
allowed orientations of the NV main axis are found and 
in general each orientation will enclose a different angle 
with the magnetic field and will be shifted by a different 
amount. 

However, if the magnetic field direction is oriented 
within the (001) plane of the diamond, always two of four 
orientations exhibit the same angle with the field and at 
special angles all four are tuned by the same amount. 

In general we therefore have to distinguish between the 
TO^^^ = ±1 states of subensemble I and II. The Hamilto- 
nian for one NV subensemble can be written as 



Hnv = ffMsB • S + DS^, + E{Sl - Si), 



(1) 



where the first term describes the Zeeman effect with 
g = 2 for an NV center and hb being the Bohr mag- 
neton. The second term denotes the zero-field splitting 
with D — 2.88 GHz and typical strain-induced E pa- 
rameters of several MHz. The presence of non-zero E 
parameters in large NV ensembles causes a mixing of the 
pure Zeeman states for low magnetic field values which 
are denoted as |0, ±). As the magnetic field amplitude is 
increased, the Eigenstates are again well approximated 
by pure Zeeman states. 

The transition frequencies lo]^^ between the eigenstates 
|0)jjj and |±)i ii are depicted in Fig.gb). 

A. Theoretical model 

Assuming that only the transition |0)j — > |— )j with 
frequency wL is in or close to resonance with the cav- 
ity mode, we can reduce the description of a single NV 
to a 2-level system. The effects caused by the presence 
of the far detuned state can be included in a con- 
stant effective frequency shift as addressed in Sec. HI A 



Htc = tUcalac + 



1 ^ 



N 



H.c. 



(2) 



where H = 1. The first two terms describe the unper- 
turbed energies of the cavity, with frequency uJc and cre- 
ation operator aj!, and of the N ensemble spins using 
the usual Pauli spin operators ([cj^jCr"] = cTiSij). Each 
spin can have a different frequency iOj which is statisti- 
cally spread around the center frequency wL. The third 
term describes the coupling to the cavity with individ- 
ual strength gj. Assuming that the ensemble spins are 
confined to a volume small compared to the wavelength, 
the ensemble will interact collectively with the mode. In 
the case of few excitations, the ensemble behaves like 
a harmonic oscillator which couples to the mode with 
the collective coupling strength ft. The collective cou- 
pling strength is given by il = yX]f l.9j l^ which for 



identical gj = (?singic = 9 gives il — gv N . To in- 
clude the probe field of the cavity we add another term 



i {rja 



?7*ace'"f*) to Eq.[2] 



Hp 

Effects of the coupling to a finite temperature bath are 
analyzed in Sec jlll B"| where we study the master equa- 
tion of the ensemble-cavity system. 

For an ensemble with identical frequencies and cou- 
pling strengths we find that we can reach the regime of 
coherent oscillations between the cavity and the ensemble 
spins if the collective coupling 17 dominates the linewidth 
of the cavity k, and the decay rate of the single spin 7hom- 
This is commonly known as strong coupling regime. In 
case of an ensemble with inhomogeneous frequency distri- 
bution it is not immediately clear under which conditions 
we can observe the avoided crossing. Here we address the 
influence of the width and form of the inhomogeneous 
frequency distribution on the avoided crossing. 



B. Experimental setup and parameters 

The experimental set-up has been described already in 
detail in [7] and is explained here only in brevity in the 
following. The heart of the experiment is a A/2 super- 
conducting coplanar waveguide (CPW) resonator with a 
center frequency of 2.7 GHz that is cooled down to 20 mK 
in a dilution refrigerator. In order to couple NV defect 
centers to the microwave field in the cavity, a (OOl)-cut 
single crystal diamond is placed in the middle of the res- 
onator, where the oscillating magnetic field exhibits an 
antinode. A meander geometry of the resonator ensures 
that a large fraction of the magnetic mode volume is cov- 
ered by spins. The high-pressure high-temperature dia- 
mond chosen in this setup contains an NV concentration 
of about 6 ppm, which corresponds to an average sepa- 
ration of about 10 nm. 
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A 2-axis Helmholtz coil configuration creates a homo- 
geneous static magnetic field oriented with an arbitrary 
direction within the (001) plane of the diamond, which 
is at the same time parallel to the resonator chip surface. 
Since perpendicular magnetic field components with re- 
spect to the resonator surface would shift the resonance 
frequency of the cavity, great care was taken during the 
alignment of the set-up. With the diamond on top of the 
chip the cavity quality factor is Q = 3200. 

In a typical experiment, we first set the magnetic field 
direction and amplitude and then measure the microwave 
transmission through the cavity with a vector network 
analyzer. 



III. CAVITY TRANSMISSION SPECTRA FOR 
INHOMOGENEOUS ENSEMBLES 

As generic experiment to test and characterize the sys- 
tem properties we analyze the weak field probe transmis- 
sion spectrum through the resonator in the weak field 
limit, where the number of excitations is negligible com- 
pared to the ensemble size. Hence, at least at low tem- 
peratures, we can largely ignore saturation effects and 
apply various simplified theoretical descriptions to ex- 
tract the central system parameters. In fact as we deal 
with more than N > 10^^ spins in a transition frequency 
range of about 10'^ Hz, we have several thousand spins 
per Hz frequency range and up to a million spins within 
the homogeneous width of at least several kHz. Hence 
theoretical modelling as a collection of effective oscilla- 
tors should provide an excellent model basis. This has to 
be taken with care at higher temperatures kgT ^ huj^, 
when we have a significant fraction of the particles ex- 
cited. 



A. Coupled collective oscillator approximation 

A photon that enters the cavity will be absorbed into a 
symmetric excitation of the ensemble spins with a weight 
given by the individual coupling strength. Subsequently 
the broad distribution of the spin frequencies induces a 
relative dephasing of this excitation so that the backcou- 
pling to the cavity mode is suppressed. As spontaneous 
decay (Ti « 44 s) is negligibly slow [7] , the decay rate of 
the polarization of the ensemble is just proportional to 
this dephasing and thus the inhomogeneous width of the 
spins. 

As the number of photons entering the ensemble is 
very small compared to iV in a first model, we simply 
approximate the ensemble as a harmonic oscillator with 
frequency uiai = wL and an effective large width 7 that 
mimics the inhomogeneity. The ensemble oscillator is 
coupled to the cavity with frequency uJc and decay rate k. 

This very simplified model already allows us to study 
the coupling dynamics of the collective energy levels of 
the broadened ensemble and the cavity mode, and in par- 



ticular the avoided crossing of the energy levels, when the 
relative energies are varied. 

To keep the model as simple as possible but still grasp- 
ing the essential properties, the presence of the other of- 
fresonant spin ensembles at w^J and uj^_^ is modelled 
by additional oscillators with frequency ujaj, (j — 2, 3, 4) 
and equal decay rate 7. 

For a small probe field injected into the cavity with 
the frequency ujp the corresponding equations in a frame 
rotating with ujp then read: 



dt 



(4) 



with probe amplitude 77, A^, = — and A^j = uJaj — 
LOp. From Eqs. [3] and |4] we can calculate the steady state 
field in the cavity (flcT^^, which can be written as 



V 

K + Ta+i{l^c- Ua) 



2 w, ' 



7/2+iA„ 



(5) 



where 



^2 I A2 



and Ua = ^ 



4 2 

9 

^.^,(7/2r + A2^. ^ ^t^(7/2r+A,^. 



2 _L A2 

(6) 



In the absence of the offresonant levels (A^2 = A3 = 
N4 — 0) we recover the situation, where for ojai = oJc 
and 5\/Ai > (7/2 — k)/2 we find two normal modes split 
by 2yJ g^Nx — (7/2 — k)2/4. The offresonant transitions 
make the situation more complex. As it can be seen 
in Eq. [5] they induce a shift Ua of the cavity frequency 
and increase the decay rate of the cavity by Fq. Both, 
shift and additional decay rate, depend on the probe fre- 
quency Wp. For the parameters in our measurement we 
find that F^ is negligible compared to n. We further ap- 
proximate \Ja by setting cjp — Wc (as the scan range of 
Wp around Wc is small compared to uj^j {j = 2, 3, 4)). We 
obtain the model of two coupled oscillators, where one of 
them has been shifted by the offresonant transitions. As 
a function of ujp and Wai (tuned by the magnetic field) , 
I (ac)st 1'^ shows an avoided crossing. Although the offreso- 
nant transitions are shifted by the magnetic field as well, 
the effect in Ua is negligibly small, so that we assume 
that Ua is constant. 

We show the measured signal at the avoided crossing 
m Fig.[l];a). By fitting |( flc) to the normal mode split- 
ting we can deduce the inhomogeneous width 7 and the 
effective coupling gV N. The Q factor of our cavity corre- 
sponds to k/ (27r) = 0.4 MHz, then 2k is the full width at 
half maximum (FWHM) of the cavity resonance. An ex- 
emplary fit result is shown in Fig. [IJc) . We have included 
a linear term in the fit to account for a background in the 
data. In Fig.[ljd) we plot |(ac)stl^ using the obtained pa- 
rameters. 
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/I [GHz] A, [MHz] /aJGHz] 

FIG. 1. (a) (Color online) Measured transmission as a function of /i = lJ_/{2-k) tuned by the magnetic field and pump 
frequency fp = u!p/{2-K). (b) Transition frequencies f]l^^ = Lj/^^/{2n) as a function of the magnetic field for a field angle of 
(fi — 22.5 . The frequency of the cavity fc = tUc/ {2tv) is denoted by the (red) dashed line, (c) Fit of | (ac)gtj^ to the transmitted 
signal close to the resonance. We obtain the decay rate 7 = 10.92 MHz. (d) (tSc)^^ as a function of fa^ and fp incorporating 
the parameters obtained from the fit. 
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FIG. 2. Simulation of the normal mode splitting at resonance 
(ojc — (7a = ijJai) for different values of 7. Finally for large val- 
ues of 7 the doublet is not resolved any more and the avoided 
crossing (sketched on the right) becomes blurred. 



From the fit we deduce a decay rate 7 = 10.92 MHz 
of the ensemble oscillator and an effective coupling of 
g\/N = 9.51 MHz. To demonstrate the effect of the en- 
semble oscillator width 7 on the normal mode splitting, 
we plot |(ac)stl^ for ujc — Ua = ujai and different values of 
7 in Fig. [2] 

Treating the ensemble as a broad harmonic oscilla- 
tor with decay rate 7 corresponds to the assumption 
that the frequency distribution of the real ensemble is a 
Lorentzian distribution, given that all spins couple with 
equal strength [13 . 



B. Polarization decay and collective coupling at 
finite temperature 

Despite cooling to very low temperatures is possible 
in the experiments, it is still important and instructive 
to study the role of thermal excitations in the system. 
In contrast to previous models based on virtually zero 
T atomic ensembles, the NV centers are in thermal con- 
tact with the chip at small but finite T. We will now 
investigate how sensitive the system reacts on thermal 
fiuct nations. 

At this point we include any shifts caused by the off- 
resonant ensembles in an effective detuning and concen- 
trate on the collective coupling between the cavity at lUc 
and the near resonant ensemble centered around ojL . For 
simplicity we assume equal coupling g for all spins to get 



H' 



TC 



ojcaiac + 



1 ^ 



9 



N 

E 

3 



H.c. 



(7) 

The effects of unequal coupling of the spins is addressed 
in [Tij. 

To include thermal excitations of the mode and the 
ensemble we have to add standard Liouvillian terms to 
the dynamics and study the corresponding master equa- 
tion of the reduced cavity ensemble system |15j . where 
the inhomogeneous width of the ensemble is still simply 
approximated by an effective dephasing term 7j, for the 
polarization. In this model decay 7hom and dephasing 
7p are described by separate quantities, which already 
should improve the model. The master equation reads 
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H 



TC " 



H, 



P 



(8) 
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where 



{crf{T))^^ can be written as [5] 



C[p\ = K. {n (T, Wc) + 1) {2acpal - a^a^p - pa\a^ 
+ Kfi {T, Wc) {2a\pac - a^alp - pUcal) 

N 



7ho 



7ho 



2 



TV 



(9) 



The first two hnes of Eq. |9] describe the couphng of the 
cavity to the bath, while the next two lines include the 
coupling of the ensemble to the bath. The number of 
thermal excitations at temperature T and frequency lu 
is denoted by n{T,uj). The term in the last line intro- 
duces nonradiative dephasing at a rate 7p of the spins 
and thereby models the inhomogeneity. 

Based on the master equation we can derive a hierar- 
chic set of equations for various system expectation val- 
ues starting with 



-- (flc) = Tr{ac--p} 



d_ 
dt 



= - (k + iAc) (ac) - igN {a^)+r] , (10) 
(-r) = - + 7hon.n (r,a.L) + 7p + IAL) {a;) 

(11) 



1.9 {<^i ac 



dt 



(af ) - -i2g ((a+a,) - (a" at)) - 7ho„, («) + 1) 



(12) 



which also includes equations for (^QcCr^), (acCTi), 
W'^j)' (4ac), {accr~), (44), {a~a~), (crfcr+), 
{^a^a^') and (a\acCr^Y In order to truncate the system 
higher order terms in the equations are expanded in a 
well defined way ^16) . then higher order cumulants are 
neglected [TTl [18] . The equations again are written in a 
frame rotating with the cavity probe frequency LOp. De- 
spite the inhomogeneous broadening, which is included 
by the nonradiative dephasing of the spins, we assume 
that all spins are equal so that we only have to include 
the equations for one spin i and pairs i, j |17j . This set of 
equations can be integrated numerically to study the dy- 
namics of the coupled system. We note that when we fix 
(erf) = — 1 we immediately arrive at the model discussed 
III Alwith 7 = 7hom 



Sec. 



■27p 



In our experiment we measure the avoided crossing for 
different temperatures of the environment and compare it 
to the results of our model. First we note that the steady 
state of the inversion as a function of the temperature 
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1 + 2n {T.ujL) 

huj- 



= tanh 



2kBT 



(13) 



For higher temperatures (cf (T))^^ is reduced and there- 
fore the effective number of NVs that take part in the 
dynamics. In the model equations this is represented 
by the last term in Eq |ll| for the polarization involving 
(flcuf ), which leads to a cutoff for the coupling at higher 
T. As a zero field approximation we thus write 



n{T) = ^WiVtanh 



hcul 
2kBT 



(14) 



where we replaced LUi by the center frequency wL . Equa- 
tion ( 14 ) should give an approximate description of the 



reduction of the Rabi splitting with increasing tempera- 
ture. 

This treatment however neglects the presence of the 
ms = +1 state which will also be populated with increas- 
ing T. Including this level and assuming n(T,u!^_^_^ « 
n (T, wL) we find the population difference between the 
ms = and ms = — 1 state to be 



A^-i(T) 
N 



NojT) 
N 



1 

1 -t- 3n (T, ujI 



This suggests 



n{T) = g 




(15) 



(16) 



to be a better description for the temperature depen- 
dence of the coupling. 

In a second step we integrate the whole hierarchic set 



of equations numerically for lOc 



and varying pump 



frequency uip. From |(ac)stl i'^) determine the Rabi 
splitting for different temperatures from which we obtain 
the collective coupling. This can be compared to our 
measurements and the approximations in Eqs. 



14 and 



m 

In Fig. [3](a) we show that the measured coupling 
strength is in disagreement with the Eq. [M] which is 
most significant for very low temperatures. Both func- 
tions Eq. [M] and Eq. [16] fail to capture the behavior for 
low temperatures. The disagreement is even more pro- 
nounced as we include the effect of the ms = +1 state. 
So far we found no definite explanation for the disagree- 
ment. However, one possible explanation would be that 
for low temperatures not all defects are "active". As 
temperature increases, more NVs become available but 
at the same time the number of NVs taking part in the 
dynamics is proportional to tanh {hu?_ / {2kBT)^ . 

The collective coupling strength determined from the 
numerical integration of the coupled equations, shown in 
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FIG. 3. (Color online) (a) Collective coupling strength ob- 
tained from our measurements for different temperatures 
(black markers) and Q{T) = ao [tanh {hujL/{2kBT))Y^^ (red 
dashed line), where ao was chosen to match the high tem- 
perature behavior. Including the ms = +1 state via Eq. [161 
gives the result denoted by the blue dotted line, (b) Collective 
coupling strength determined from the results of the numer- 
ical integration of the coupled equations with g^/N — ao (ao 
determined from (a), black markers) and n(r) (red dashed 
line). 



Fig. ^h), exactly follows Eq. 
assumption n{T) oc (crf(r)),t 
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This shows that the 
is reasonable. However, 
the almost constant value of the coupling strength il{T) 
found in the experimental data for very low T cannot be 
explained in the above theoretical model. As possible ex- 
planations one might think of a reduced thermal excita- 
tion probability due to spin-spin coupling or an effective 
reduction of active NV centers close to zero temperature. 
Interestingly the same behavior is also found in measure- 
ments of dispersive shift of the cavity mode as a function 
of temperature by the offresonant spin ensemble at zero 
magnetic field. 



C. Detailed modeling and reconstruction of 
inhomogeneous distributions 

The simplified model descriptions discussed above 
in Sees. Ill A| and |III B| provided for an analytically 
tractable and qualitatively correct description of the ef- 
fect of an inhomogeneous broadening of the ensemble. 
This also allows to get a fairly good estimate for the 



total width of the frequency distribution of the ensem- 
ble. However, such an effective width model inherently is 
connected to the assumption of a Lorentzian shape of the 
ensemble frequency distribution. In actual crystals such 
an assumption is not obvious and other distributions of 
local field variations and strain distributions are possible 
as well. 

To obtain more accurate information about the dis- 
tribution we will now use an improved model based on 
the resolvent formalism to treat the coupling between a 
central oscillator (the mode) and the spin degrees of free- 
dom [131 [in]. Here each frequency class of spins is treated 
individually. For low temperatures virtually all spins are 
in the ms — state, i.e. the lower state of our effective 
two-level system and their excitation properties can be 
approximated by a frequency distributed set of oscillators 
(Holstein-Primakoff approximation) . 

We define creation and annihilation operators for the 
corresponding ensemble oscillators representing a sub- 
class of two-level systems with equal frequency via 



-1 



2 a] a j 



and 



1 



a] a J « a] 



(17) 

The approximation in Eq. [17] is justified as long as the 
number of excitations in each ensemble is much smaller 
than the number of spins in this energy region. In our 
experimental setup this is very well justified and we thus 
obtain the unperturbed part of the Hamiltonian as 



N 



Hg = LJcdlcLc + ^ 



and the interaction term 

N 

V = ^ (gja]ac + 9*ajal] , 



(18) 



(19) 



which constitute H' = Hq+V'. In this section we ac- 
count for the decay of cavity excitations and the spon- 
taneous decay of the spins by introducing nonzero imag- 
inary parts of the corresponding transition frequencies 
ImUc = —K and Imujj — — ^7hom and the resolvent of 
the Hamiltonian H' is defined as G (z) = 1/ (z — H'). 

Let us consider the state \(pc) = |lc, 0, . . . , 0) where we 
have one photon in the cavity and no excitation in the 
ensemble. 

The matrix element of the resolvent Gcc (z) — 
{ifcl G (z) \ipc) can be written as 

where we define the matrix element of the level shift op- 
erator 



Z H/i 



y y K^^v:,^i^vL + • ■ ■ , (21) 
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where V^i = {ipk\V \(pi) and = Hp |(^c)- States 
\ipi) with i 7^ c are states where the excitation is absorbed 
in the ensemble spin i. We note that — and that 
V^'j = for i, j ^ c since our Hamiltonian does not include 
spin-spin interaction. Only the second term in Eq. 21 
remains and we can write 



Rcc (Z) = J2 



z-E,, 



(22) 



Introducing the coupling density profile p (lo) 
J2j \9j\ S {uj — Reujj) as it is done in [T3j leads to 



Rcc (z) 



p{uj) 



-duj 



>7ho 



(23) 



Approaching the branch cut at w — h^hom we write 



lim RcciuJ- ^7hom ± i?? ) = Rfc i ^ - ^7hom 

T]^o+ \ 2 J \ 2 



(24) 



Using lim,,^o+ = '^x'^^'^^ (^)' where V denotes the 
Cauchy principal value, we find 



V 



Pico') 



ydw' T iTTp (w) . (25) 



Experimentally we probe the transmission of a weak 
probe signal amplitude through the cavity as a function 
of frequency. The position and shape of the weak field 
transmission resonances can be determined from the com- 
plex poles of G+ (w) = 1/ (cl> — Wc — Rtc (^)); which con- 
tains the spin energy distribution on the right hand side. 
We can therefore extract information about the coupling 
density p{uj) = — ^Im_R+ (w — |7hom) by carefully an- 
alyzing the measured transmission spectrum. For small 
7hom the reconstruction simplifies to: 



p{uj) 



1 



1 aRei?+(w) 



27r 



duj 



-7ho 



(26) 



where Imi?+ (u) — 57hom) is expanded in a Taylor series 
around 7hoin = 0. We can therefore directly use the fre- 
quency distribution of the transmitted signal |(ac)g^|^ via 



[uj-ujc- Rci?+ (cj))^ + + Imi?+ {u:)f 

(27) 

to determine i?+(a;). Let us point out here that 
Eq. [27] exhibits a Lorentzian shape as a function of 
with peak position w — Re i?+ [uj) and peak width of 
2 (k -I- Imi?+ (w)). After extraction of the relevant pa- 
rameters which determine i?^ (oj) from the measured 
spectra, we can simply use Eq. [26] to find the coupling 
density p{uj). Assuming that all spins are coupled with 
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FIG. 4. (Color online) Coupling density of the ensemble de- 
termined from an exemplary transmission measurement at 
T = 20 mK. The errorbars are calculated from the uncer- 
tainties in (ij) which we obtain from fittin g Eq. 
data. We fit the q-Gaussian (red) in Eq. 



28 



27 



to the 



to the coupling 
density to determine the width and the behavior of the tails 
of the distribution. 



equal strength, one finds that p (w) = g^NAi {uj) where 
M (uj) is the frequency distribution of the spins. 

The shape of the coupling density, i.e. the frequency 
distribution of the spins, plays an important role in the 
cavity ensemble interaction [T^l [13] . If the coupling den- 
sity falls off sufficiently fast with distance from the center, 
the width of the Rabi peaks will decrease with increasing 
collective coupling strength g\/N . For the spin frequency 
distribution, the limiting case is the Lorentzian coupling 
density profile, for which the width of the Rabi peaks 
is independent of g\fN [12]. Any distribution falling of 
faster than will provide a decrease of the width of 
the Rabi peaks. Moreover, knowing the coupling den- 
sity gives us the opportunity to study the transmission 
through the cavity for different parameter ranges via 
Eq.[2§ 

To determine p {lo) the raw data have to be rearranged, 
since in the experiment we cannot simply vary ujc but 
we shift the center frequency of the spins by a magnetic 
field. We therefore shift each scan by — [B) to ob- 
tain fixed ensemble frequencies and tuning of the cavity 
frequency. As the only significant quantity is the detun- 
ing between the cavity and the ensemble this does not 
change the dynamics. For fixed ui we fit a Lorentzian to 
|G+ (w, Wc)|^ to determine i?+ (w) in order to calculate 
p{uj). We plot an exemplary result for p (w) in Fig. ]4] 
To determine the behavior of the tails of the distribution 
we fit the function 



[lo - cjo) 



(28) 



to the data. L (w) is related to the q-Gaussian, a Tsallis 
distribution. The dimensionless parameter 1 < g < 3 
determines how fast the tails of the distribution fall 
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FIG. 5. (Color online) In (a) we show the coupling density 
for different temperatures. The g-parameter and width 7, 
of the coupling density do not change significantly with in- 
creasing temperature. The coupling g\fN = (J p (oj) do;)^''^ 
is reduced, as can be seen in (b). The red line shows 

f2(r) = ao [tanh (ftajl_/(2fcsT))] ^ with ao chosen to match 
the high temperature limit of the data (black markers). To 
compare ^{T) with our previous results determined directly 
from the splitting we plot the data already shown in Fig. |3] 
using grey markers. 



off, while a is related to the width. The actual width 
(FWHM) is given by 7,, = '^^ °'%Zl^ ■ For g ^ 1 we 
recover a Gaussian distribution, while for q — 2 we find a 
Lorentzian distribution. The wings fall of as 

From the fit in Fig. [4] we find values for q = 1.389 ± 
0.007 and 7q/(27r) = (12.54 ± 0.10) MHz. The tempera- 
ture during the measurement was 20 mK. 

The same analysis is performed for data measured for 
T 20 - 1000 mK, see Figjsja). The fitted curves are 
not shown. For the resulting coupling densities we find 
almost no change in the width or the q-parameter. How- 
ever, with increasing temperature J p (w) do; = g^N is 
reduced, as can be seen in Fig. [Sjb). The resulting cou- 
pling strength is in agreement with the results obtained 



from the analysis in Sec HIB It also reproduces the un- 
expected behavior for small values of T. 

To study the behavior of an ensemble following a q- 
Gaussian distribution with the q-parameter we deter- 
mined, we show |G'cc(i^)| for w^, ~ ujq as & function of 
the collective coupling strength g^fN in Fig. gb). The 
transmission for an ensemble with Gaussian distribution 




d) qVN — 4MHz 









-30 -20 -10 

/|MHz] 



■■) gVN = SMIIz 



-20 



-10 10 

/[MHz] 



20 30 





- 16MHz 




1 ; 
1 : 

1 ! 

1 ■■ 


a 











30 -20 



-10 

/[MHz] 



20 30 



FIG. 6. (Color online) We plot \Gcc (w)|^ as a function 
of the collective coupling strength gy/N and for three dif- 
ferent values of the q-Gaussian parameter q (grayscale fig- 
ures on the left are plotted on a logarithmic scale). For 
g\fN = 4, 8 and 16 MHz, marked by d), e) and f) we show 
the transmission for all three values of q together. The trans- 
mission for the Gaussian coupling density with g — >■ 1 (dot- 
ted black line), for the coupling density of our ensemble with 
q = 1.39 (dashed red line) and for the Lorentzian coupling 
density with q — >■ 2 (solid black line) are shown on the 
right. The parameters were chosen to 7q/(27r) = 10 MHz, 
K/(27r) = 0.4 MHz and 7hom/(27r) = 1 Hz. 



(g — !■ 1) is shown in Fig. J6^ a) and for a Lorentzian dis- 
tribution {q — > 2) in Fig. ^c). The width parameter a is 
chosen accordingly to ensure that 7^ is the same for all 
three cases. In Figs.[6] d)-f) the transmission for the three 
ensemble types is shown for g\/~N = 4, 8 and 16 MHz, 
respectively. For the ensemble with Lorentzian distri- 
bution (solid black line) we see that the width of the 
Rabi resonances remains constant with increasing col- 
lective coupling. For the Gaussian distribution (dot- 
ted black line) and the intermediate distribution with 
q = 1.39 (dashed red line) we find a decrease in the 
peak width for increasing collective coupling. In Fig. [7] 
we plot the real and imaginary part of one of the com- 
plex poles of G+j (w) , determining position and width of 
one of the Rabi peaks. As we focus on the resonant 
case the spectrum is symmetric. We chose 7hom = 0, 
K/(27r) = 0.4 MHz, 75/(27r) = 10 MHz and q = 1.39. 
This again shows the decreasing width of the Rabi Peaks 
as gyfN increases. 

We therefore assume that for our ensemble it is possible 
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FIG. 7. (a) Real part and (b) modulus of the imaginary part 
of one of the complex poles of (cj) in the regime where 
gVN/(2-n-) > 7g/(27r) = 10 MHz. The collective couphng 
strength g\fN is shown as dotted line in (a). With increasing 
g^/N the modulus of the imaginary part, proportional to the 
width of the Rabi peaks, is reduced. 



to increase the lifetime of the collective states by increas- 
ing the collective coupling. This could be achieved by a 
further decrease of the mode volume or an increase of the 
NV density in the sample. 

In Fig. [8] we show the behavior of the transmission 
with increasing 7^. For the Lorentzian coupling density 
the splitting of the resonance peaks is always reduced if 
7g > 0. For coupling densities falling off faster than 
as it is the case in (a) and (b), the splitting of the reso- 
nance peaks is even slightly increased for 7g > 0, until the 
peaks finally merge. Hence a finite but small enough in- 
homogeneous width might even mimic somewhat higher 
active spin numbers. 
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FIG. 8. \Gcc (i^)!^ on resonance as a function of the inhomo- 
geneous width 7^ for an ensemble with (a) Gaussian coupling 
density, (b) q-Gaussian (g = 1.39) coupling density and (c) 
Lorentzian coupling density. To keep the transmission visible 
for large values of 7, the range of the color coding is limited to 
[0,0.1]. The parameters were chosen to (;\/iV/(27r) = 10 MHz, 
K/(27r) = 0.4 MHz and 7hom/(27r) = 1 Hz. 



IV. A TRANSMISSION LINE MICRO-MASER 
WITH AN INHOMOGENEOUS NV ENSEMBLE 



As bottom line we see that common coupling to a cav- 
ity mode can suppress dephasing of the polarization, if 
the effective Rabi frequency is larger than the inhomoge- 
neous width. This could be interpreted as the effect that 
the exchange of the excitation between the ensemble and 
the cavity is so fast, that there is no time for decoherence 
in the ensemble. 



A recent proposal to construct a laser operating on an 
ultra-narrow atomic clock transition predicted very nar- 
row optical emission above threshold [T7]. Similar ideas, 
employing the collective coupling between a cold atomic 
ensemble and a microwave cavity, have been proposed 
to construct stable stripline oscillators in the microwave 
regime [TS] . 
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Here we study the prospects of implementing such an 
oscillator by coupling a diamond to the CPW resonator. 
At first sight in view of the MHz scale inhomogeneous 
broadening, one would expect fast dephasing. However, 
as we have seen above, for strong enough coupling one ob- 
serves a continuous rephasing of the polarization inducing 
a long lived polarization and coherent Rabi oscillations. 
Thus one could expect narrow microwave emission nev- 
ertheless. 

Let us consider the case of the cavity mode tuned to 
resonance with the spin transition |0)j — >■ |— )j, which 
is partially inverted by an external incoherent pump. 
Such a pump could in principle be facilitated by opti- 
cal pumping and it can be consistently modeled by a re- 
versed spontaneous decay. Alternatively one could think 
of pulsed inversion by tailored microwave pulses or a time 
switching of the magnetic bias field. 



mg 



be modelled by adding 



the terms 
where w 



can 

-f Ef=i (a-<p + pa-a; -2a+pa-), 
denotes the pump rate, to the Liouvillian in Eq. [9) 

For the explicit calculations, here we use the effec- 
tive linewidth model, as outlined in Sec. |IH B| without 
any coherent pump. Hence the total phase symme- 
try of the system is not broken and we assume (oc) = 
(aj!^ — {cr^) = 0. Starting from the master equation in 
Eq. 8 we derive four coupled equations for (ajac), 
(flcOY ) ^-iid (^(T^ajY We used the cumulant expansion 

{a\ac<Ji) — i^c'^'^'^i) cuia ~^ i'^c^c) (cf) which takes this 
simple form because of the total phase invariance. As- 
suming that the higher order cumulant (^alaccrl 



) can 

/ cum 



Mathematically such 



be neglected, we arrive at a closed set of four equations 
that can be solved analytically. To study the spectrum 
of the emitted light we calculate the two-time correlation 
function (^a^ (r) Qc (0)) via the quantum regression theo- 
rem. We switch to a frame rotating with ujc and define 



incoherent 



pump- 



A = 



We obtain 



igN 

7hom"- (r, wL) 



■Ip 



iA 



(at (r)a,(0)) 
(^+W«c(0)) 



(29) 



r 



From Eq. 29 we can calculate the spectrum via Laplace 
transformation [20]. We show the linewidth of the ob- 
tained spectrum in the resonant case (A = 0) for different 
values of the pump w and the inhomogeneous width 7p , 
see Fig. [9j The minimum linewidth is A/ « / n. For 
small 7p the region where we find narrow linewidth emis- 
sion is characterized by 7hom < w < 2g'^N / k. The critical 
width of the inhomogeneity is given by 7p,crit = g'^N/n. 
Hence once we achieve enough pumping and coupling 
strength the system could provide for an extremely stable 
microwave oscillator. 



CONCLUSIONS 



tribution. Interestingly for our dense NV ensemble, the 
frequency distribution of the spins can be well described 
as a q-Gaussian with q — 1.39, showing that the wings of 
the distribution fall off faster than l/w^. The tempera- 
ture dependence of the effective available spin number fits 
quite well with expectations, except for an unexpected 
decrease for very low temperature. In summary such an 
NV ensemble cavity QED system exhibits a prolonged 
lifetime of the eigenstates of the coupled cavity-ensemble 
system JJ3J and has great potential as quantum interface 
between superconducting and optical qbits. The long ef- 
fective Ti time could also be the basis of building a com- 
pact ultrastable microwave oscillator if the strong cou- 
pling overcomes the dephasing from the inhomogeneous 
broadening. 



We showed theoretically and experimentally that an 
ensemble of spins with an inhomogeneous frequency dis- 
tribution coupled to a cavity mode can exhibit strong 
coupling, where the coherent energy exchange between 
mode and ensemble dominates cavity decay and polar- 
ization dephasing. A detailed theoretical modeling con- 
necting probe transmission and frequency distribution al- 
lows to extract not only the effective coupling strength 
and particle number but also the detailed frequency dis- 
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